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§ ■ 1 Introduction 

o 

Recent measurements of the CP- violating asymmetry parameter sin2/3 by the BABAR |J 
and BELLE || collaborations established CP violation in the B° system. These 
; measurements, as well as updated preliminary results [§, ||, are consistent with the 

Standard Model expectation based on measurements and theoretical estimates of the 
elements of the Cabibbo-Kobayashi-Maskawa || (CKM) quark-mixing matrix. 

The study of B decays to charmless hadronic two-body final states will yield im- 
portant information about the remaining angles (a and 7) of the Unitarity Triangle. 
In the Standard Model, the time-dependent CP-violating asymmetry in the decay 
■ B° — > 7r + 7r~ is related to the angle a, and ratios of branching fractions for various tttt 

and Ktt decay modes are sensitive to the angle 7. We previously reported measure- 
ments of branching fractions |J and CP- violating asymmetries |7| in B° — > 7r + 7r~, 
K + n~, and K + K~ decays. (Unless explicitly stated, charge conjugate decay modes 
are assumed throughout this paper.) In this paper, we present a preliminary update 
of these results using a sample of 60 million BB pairs. 

We reconstruct a sample of B mesons (B rec ) decaying to the h + h'~ final state, 
where h and h' refer to it or K, and examine the remaining charged particles in each 
event to "tag" the flavor of the other B meson (B tag ). The decay rate distribution 
/+ (/-) when h + h'~ = 7r + 7r~ and B tag = B° (B°) is given by 

e -|At|/r 

/±(At) = — [1 ± S m sm(Am d At) t C wn cos(Am d At)] , (1) 

4r 
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where r is the mean B° lifetime, Am^ is the eigenstate mass difference, and At = 
^rec — ^tag is the time between the B rec and B tag decays. The CP- violating parameters 
(Stttt cLlld C-ft-ft are defined in terms of a complex parameter A as 

^' 7r7r = 7772 an< ^ = ~ iT72" (^) 

1 -\- A 1 -4- A 

If the decay proceeds purely through the 6 — > wW 7- tree process, then A is given in 
terms of the CKM elements Vy by 

In this case, C nn = and S n7T = sin2o;, where a = arg [— VtaKb/KidKib]- m general, 
the b ^ dg penguin amplitude modifies both the magnitude and phase of A, so that 



C wn 7^ and S W7r = J 1 — sin 2a eff , where a e ff depends on the magnitudes and rel- 
ative strong and weak phases of the tree and penguin amplitudes. Several approaches 
have been proposed to obtain information on a in the presence of penguins [Efl. 



2 The BABAR detector and dataset 

The data sample used in this analysis consists of 55.6 fb _1 , corresponding to 60.2±0.7 
million BB pairs, collected on the T(4S) resonance with the BABAR detector at 
the SLAC PEP-II storage ring between October 1999 and December 2001. Equal 
branching fractions for T(4S) — > B°B° and B + B~ are assumed. 

A detailed description of the BABAR detector is presented in |5). Charged particle 
(track) momenta are measured in a tracking system consisting of a 5-layer double- 
sided silicon vertex tracker (SVT) and a 40-layer drift chamber (DCH) filled with a gas 
mixture of helium and isobutane. The SVT and DCH operate within a 1.5 T supercon- 
ducting solenoidal magnet. Photons are detected in an electromagnetic calorimeter 
(EMC) consisting of 6580 CsI(Tl) crystals arranged in barrel and forward endcap sub- 
detectors. The flux return for the solenoid is composed of multiple layers of iron and 
resistive plate chambers for the identification of muons and long-lived neutral hadrons. 
Tracks from the B rec decay are identified as pions or kaons by the Cherenkov angle 
9 C measured with a detector of internally reflected Cherenkov light (DIRC). 



3 Analysis method 

Event selection is identical to that described in fffl. Candidate B rec decays are re- 
constructed from pairs of oppositely-charged tracks forming a good quality vertex, 
where the B rec four-vector is calculated assuming the pion mass for both tracks. We 
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require each track to have an associated 6 C measurement with a minimum of six 
Cherenkov photons above background, where the average is approximately 30 for 
both pions and kaons. Protons are rejected based on 9 C and electrons are rejected 
based on dE/dx measurements in the tracking system, shower shape in the EMC, 
and the ratio of shower energy and track momentum. Background from the reaction 
e + e~ — > qq (q = u, d, s, c) is suppressed by removing jet-like events from the sample: 
we define the center-of-mass (CM) angle 8s between the sphericity axes of the B can- 
didate and the remaining tracks and photons in the event, and require [cos^sl < 0.8, 
which removes 83% of the background. The total efficiency for signal events of the 
above selection is approximately 38%. 

Signal decays are identified kinematically using two variables. We define a beam- 
energy substituted mass m ES = \J E^ — p^, where the B candidate energy is defined 
as E h = (s/2 + Pi • Pb)/ Ei, y/s and E^ are the total energies of the e + e~ system in the 
CM and laboratory frames, respectively, and pj and p# are the momentum vectors 
in the laboratory frame of the e + e~ system and the B mc candidate, respectively. 
Signal events are Gaussian distributed in mES with a mean near the B mass and a 
resolution of 2.6MeV/c 2 , dominated by the beam energy spread. The background 
shape is parameterized by a threshold function with a fixed endpoint given by 
the average beam energy. 

We define a second kinematic variable AE as the difference between the energy 
of the B rec candidate in the CM frame and y/s/2. Signal 7T7t decays are Gaussian 
distributed with a mean value near zero. For decays with one (two) kaons, the dis- 
tribution is shifted relative to tttt on average by — 45MeV (— 91MeV), respectively, 
where the exact separation depends on the laboratory momentum of the kaon(s). The 
resolution on AE is approximately 26 MeV and is validated in large samples of fully 
reconstructed B decays. The background is parameterized by a quadratic function. 

Candidate h + h'~ pairs selected in the region 5.2 < m-Es < 5.3 GeV/c 2 and \AE\ < 
0.15 GeV are used to extract yields and CP- violating asymmetries with an unbinned 
maximum likelihood fit. The total number of events in the fit region satisfying all of 
the above criteria is 17585. 

To determine the flavor of the B tag meson we use the same 5-tagging algorithm 
used in the BABAR sin2/3 analysis The algorithm relies on the correlation between 
the flavor of the b quark and the charge of the remaining tracks in the event after 
removal of the B rec candidate. We define five mutually exclusive tagging categories: 
Lepton, Kaon, NT1, NT2, and Untagged. Lepton tags rely on primary electrons and 
muons from semileptonic B decays, while Kaon tags exploit the correlation in the 
process b — > c — > s between the net kaon charge and the charge of the b quark. 
The NT1 (more certain tags) and NT2 (less certain tags) categories are derived from 
a neural network that is sensitive to charge correlations between the parent B and 
unidentified leptons and kaons, soft pions, or the charge and momentum of the track 
with the highest CM momentum. The addition of Untagged events improves the 
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Category 


e(%) 


D{%) 


AD (%) 


Q(%) 


Lepton 


11.1 ±0.2 


82.8 ± 1.8 


-1.2 ±3.0 


7.6 ±0.4 


Kaon 


34.7 ±0.4 


63.8 ± 1.4 


1.8 ±2.1 


14.1 ±0.6 


NT1 


7.6 ±0.2 


56.0 ±3.0 


-2.7 ±4.7 


2.4 ±0.3 


NT2 


14.0 ±0.3 


25.4 ±2.6 


9.4 ±3.8 


0.9 ±0.2 


Untagged 


32.6 ±0.5 








Total Q 


25.0 ±0.8 



Table 1: Tagging efficiency e, average dilution D = 1/2 (D B o ± Dg ), dilution differ- 
ence AD = Dgo — -Dgo, and effective tagging efficiency Q for signal events in each 
tagging category. The values are measured in the B^v sample. 



signal yield estimates and provides a larger sample for determining background shape 
parameters directly in the maximum likelihood fit. 

The quality of tagging is expressed in terms of the effective efficiency Q = J2 C e cD^, 
where e c is the fraction of events tagged in category c and the dilution D c = 1 — 2w c 
is related to the mistag fraction w c . Table [I] summarizes the tagging performance 
in a data sample £>fl av of fully reconstructed neutral B decays into D^*>~h + (h + = 
n + ,p + ,a,i) and J/ifiK* (K*° — > K + tc~) flavor eigenstates. We use the same tagging 
efficiencies and dilutions for signal nn, Kir, and KK decays. Separate background 
efficiencies for each species are determined simultaneously with S n7T and C n7T in the 
maximum likelihood fit. 

The time difference At is obtained from the measured distance between the z 
positions of the B rec and B t3jg decay vertices and the known boost of the e + e~ system. 
The z position of the B tag vertex is determined with an iterative procedure that 
removes tracks with a large contribution to the total \ 2 - An additional constraint 
is constructed from the three-momentum and vertex position of the B rec candidate, 
and the average e + e~ interaction point and boost. For 99.5% of candidates with a 
reconstructed vertex, the rms Az resolution is 180 /im (1.1 ps). We require \At\ < 
20 ps and a At < 2.5 ps, where o"At is the error on At. The resolution function for 
signal candidates is a sum of three Gaussians, identical to the one described in [[|, 
with parameters determined from a fit to the -Bfl av sample (including events in all five 
tagging categories). The background At distribution is parameterized as the sum of 
an exponential convolved with a Gaussian, and two additional Gaussians to account 
for tails. A common parameterization is used for all tagging categories, and the 
parameters are determined simultaneously with the CP parameters in the maximum 
likelihood fit. We find that 86% of background events are described by an effective 
lifetime of about 0.6 ps, while tails are described by 12 (2)% of events with a resolution 
of approximately 2 (8) ps. 

Discrimination of signal from light-quark background is enhanced by the use of 
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a Fisher discriminant T f|. The discriminating variables are constructed from the 
scalar sum of the CM momenta of all tracks and photons (excluding tracks from 
the B rec candidate) entering nine two-sided 10-degree concentric cones centered on 
the thrust axis of the B Tec candidate. The distribution of T for signal events is 
parameterized as a single Gaussian, with parameters determined from Monte Carlo 
simulated decays and validated with B~ — ► D n~ decays reconstructed in data. The 
background shape is parameterized as the sum of two Gaussians, with parameters 
determined directly in the maximum likelihood fit. 

Identification of h + h'~ tracks as pions or kaons is accomplished with the Cherenkov 
angle measurement from the DIRC. We construct Gaussian probability density func- 
tions (PDFs) from the difference between measured and expected values of 9 C for the 
pion or kaon hypothesis, normalized by the resolution. The DIRC performance is 
parameterized using a sample of D* + — > D°tt + , D° — > K~tt + decays, reconstructed 
in data. Figure |l] shows the typical separation between pions and kaons, which varies 
from 8cr at momenta of 2 GeV/c to 2.5a at 4 GeV/c, where a is the average resolution 
of6 r . 



M 7.5 



3.5 4 
Momentum (GeV/c) 



Figure 1: Variation of the separation between the kaon and pion Cherenkov angles 
with momentum, as obtained from a control sample of D* + — > D°it + , D° —* K~n + 
decays reconstructed in data. 



4 Results 

We use unbinned extended maximum likelihood fits to extract yields and CP pa- 
rameters from the B rec sample. The likelihood for candidate j tagged in category 
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Mode 


€T (%) 


N s (Events) 


B (10- 6 ) 


7T + 7T~ 


38.5 ±0.7 


104+16+7 


5.4 ±0.7 ±0.4 


K+TT- 


37.6 ±0.7 


403 ±24 ±15 


17.8 ±1.1 ±0.8 


K+K- 


36.7 ±0.7 


< 16 (90% C.L.) 


< 1.1 (90% C.L.) 



Table 2: Summary of results for total detection efficiencies ey, fitted signal yields N s 
and measured branching fractions B. 



c is obtained by summing the product of event yield n i} tagging efficiency e iiC , and 
probability 7\ c over the eight possible signal and background hypotheses % (referring 
to mi, K + tt~, K~tt + , and KK decays), 



exp 



E 



n 



3 L J 



(4) 



For the X+7r ± components, the yield is parameterized as rij = Nk^{1 ± Ak-k) /2, 
where iV^ = iV^-^ ± A^ x+7r - and A Kn = {N K -^+ - N K + v -)/{N K -^ + N K+7T -). 
The probabilities 7-^ jC are evaluated as the product of PDFs for each of the independent 
variables Xj = {rriEs, AE, T, + , 9~ , At}, where 6* + and 6~ are the Cherenkov angles 
for the positively and negatively charged tracks. We use identical PDFs for # + and 
6~ . The total likelihood £ is the product of likelihoods for each tagging category and 
the free parameters are determined by minimizing the quantity — ln£. 



4.1 Time-independent fit 

In order to minimize the systematic error on the branching fraction measurements, 
we perform an initial fit without tagging or At information. A total of 16 parame- 
ters are varied in the fit, including signal and background yields (6 parameters) and 
asymmetries (2), and parameters for the background shapes in itles (1), AE (2), and 
^(5). 

Table ^ summarizes results for total efficiencies, signal yields and branching frac- 
tions. The upper limit on the signal yield for B° — > K + K~ is given by the value of n° 
for which J n C max dn/ / °° C max dn = 0.90, where jC max is the likelihood as a function 
of n, maximized with respect to the remaining fit parameters. The branching frac- 
tion upper limit is calculated by increasing the signal yield upper limit and reducing 
the efficiency by their respective systematic errors. The fit result for the Ktt charge 
asymmetry Akw is 

A K7T = -0.05 ±0.06 ±0.01, 90% C.L. - 0.14 < ^ < 0.05. (5) 

The statistical and systematic errors on Ak-k are added in quadrature when calculat- 
ing the 90% confidence level (C.L.). 
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The dominant systematic error on the branching fraction measurements is due to 
uncertainty in the shape of the 9 C PDF, while the dominant error on Ak-k is due to 
possible charge bias in track and 9 C reconstruction. All measurements are consistent 
with our previous results reported in ||. 

Figure || shows distributions of tties and AE after a cut on likelihood ratios. We 




5.2 5.225 5.25 5.275 5.3 -0.1 0.1 



Figure 2: Distributions of ttj-es (left) and AE (right) for events enhanced in signal 
7T7T (top) and Kir (bottom) decays based on the likelihood ratio selection described in 
the text. Solid curves represent projections of the maximum likelihood fit result after 
accounting for the efficiency of the additional selection, while dashed curves represent 
qq and <-> Kit cross-feed background. 



define TZ sig = J2 s n s V s / 'Ei^i^t and K k = n k V k j E S "A where £ s (E») indicates a 
sum over signal (all) hypotheses, and V k indicates the probability for signal hypoth- 
esis k. The probabilities include the PDFs for 9 C , J 7 , and ities (AE) when plotting 
AE (rn.Es)- The selection is defined by optimizing the signal significance with respect 
to 7l S i g and lZ k . The solid curve in each plot represents the fit projection after cor- 
recting for the efficiency of the additional selection (approximately 67% for im and 
88% for Kir). 
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4.2 Time-dependent fit 



The time-dependent CP asymmetries and C n7T are determined from a second fit 
including tagging and At information, with the yields and Ak-k fixed to the results 
of the first fit. The At PDF for signal 7r + 7r~ decays is given by Eq. ([[]), modified to 
include the dilution and dilution difference for each tagging category, and convolved 
with the signal resolution function. The At PDF for signal Kir events takes into 
account B°-B° mixing, depending on the charge of the kaon and the flavor of B ta . g . We 
parameterize the At distribution in B° — > K + K~ decays as an exponential convolved 
with the resolution function. 

A total of 34 parameters are varied in the fit, including the values of S„ n and C nn , 
separate background tagging efficiencies for inr, Kir, and KK events (12), parameters 
for the background At resolution function (8), and parameters for the background 
shapes in m ES (5), AE (2), and T (5). The signal tagging efficiencies and dilutions 
are fixed to the values in Table [I], while r and Am^ are fixed to their PDG values [E^j. 



For each parameter, we also calculate the 90% CL. interval taking into account the 
systematic error. The fit yields 

S wn = -0.01 ± 0.37 ± 0.07, 90% CL. - 0.66 < S nw < 0.62, 

C W7T = -0.02 ± 0.29 ± 0.07, 90% CL. - 0.54 < C W7V < 0.48, (6) 

and the correlation between S nn and C n7r is —13%. 

Systematic uncertainties on S nn and C nn are dominated by the uncertainty on the 
shape of the 8 C PDF. Since we measure asymmetries near zero, multiplicative sys- 
tematic errors have also been evaluated (0.05). We sum in quadrature multiplicative 
errors, evaluated at one standard deviation, with the additive systematic uncertain- 
ties. 

To validate the analysis technique, we measure r and Am^ in the B Tec sample 
and find r = (1.66 ± 0.09) ps and Arrid = (0.517 ± 0.062)ft ps _1 . Figure |3] shows 
the asymmetry A mix = (N unmixcd - AU xed )/(AU mixed + AUixcd) in a sample of events 
enhanced in B —>■ Kir decays. The curve shows the expected oscillation given the 
value of Arrtd measured in the full sample. 

For tagged events enhanced in signal im decays, Figure [| shows the At distribu- 
tions and the asymmetry A^^At) = [N B o(At)-N m (At)]/[N B o(At)+N go (At)}. The 
selection procedure is the same as for Figure 0, with the likelihoods defined including 
the PDFs for 9 B , m ES , and AE. 



5 Summary 

In summary, we have presented updated preliminary measurements of branching frac- 
tions and CP- violating asymmetries in B° — > 7r + 7r~, K + tt~ , and K + K~ decays. All 
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Figure 3: The asymmetry A m i x between unmixed and mixed events in a sample 
enhanced in decays. The curve indicates the expected oscillation corresponding to 
Am d = 0.517 hps" 1 . The dilution from qq events is evident in the reduced amplitude 
near \At\ =0. 



results are consistent with our previous measurements. No evidence for CP violation 
is observed and our measurement of Ak-k disfavors theoretical models that predict a 



large asymmetry [13], [14 



We are grateful for the extraordinary contributions of our PEP-II colleagues in 
achieving the excellent luminosity and machine conditions that have made this work 
possible. 
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